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ABSTRACT 

Here we present a self-consistent, bimodal stationary solution for spherically 
symmetric flows driven by young massive stellar clusters with a central super- 
massive black hole (SMBH). We demonstrate that the hydrodynamic regime of 
the flow depends on the location of the cluster in the 3D (star cluster mechani- 
cal luminosity - BH mass - star cluster radius) parameter space. We show that 
a threshold mechanical luminosity (L cr i t ) separates clusters which evolve in the 
BH dominated regime frome those whose internal structure is strongly affected 
by the radiative cooling. In the first case (below the threshold energy) gravity of 
the BH separates the flow into two distinct zones: the inner accretion zone and 
the outer zone where the star cluster wind is formed. In the second case (above 
the critical luminosity), catastrophic cooling sets in inside the cluster. In this 
case the injected plasma becomes thermally unstable that inhibits a complete 
stationary solution. 

We compared the calculated accretion rates and the BH luminosities with 
those predicted by the classic Bondi accretion theory and found that Bondi's 
theory is in good agreement with our results in the case of low mass clusters. 
However, it substantially underestimates the accretion rates and BH luminosities 
if the star cluster mechanical luminosity, L$c > 0.1L crit . 

Subject headings: accretion — galaxies: active — galaxies: starburst — galaxies: 
star clusters — hydrodynamics 



1. Introduction 



Intensive studies of active galactic nuclei (AGNs) in the optical, infrared (IR) and X- 
ray regimes during the last decade, have unveiled the presence of massive starbursts around 
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the central supermassive black hole (BH) in a number of Seyfert galaxies. For example, 
Rodriguez Espinosa et al. (1987) found that the 25, 60 and 100 /im fluxes of classical 
optically selected Seyfert galaxies are not correlated with their ultraviolet (UV) and optical 
continuum and that far-IR colors of the selected galaxies are indistinguishable from those of 
starburst galaxies. They suggested then that high far-IR luminosities associated with many 
Seyfert galaxies indicate an intrinsic link between the circumnuclear star formation and the 
AGN activity. Baum et al. (1993) revealed a kiloparsec-scale, diffuse radio emission in 12 of 
the 13 observed Seyfert galaxies. They found that the intensity of the diffuse radio emission 
correlates with the far-IR luminosity of the host galaxy, suggesting that this emission is 
generated in the galactic superwind plasma driven by the nuclear starburst along the minor 
axis of the galaxy, and claimed that circumnuclear starbursts and starburst driven outflows 
may be intrinsic to many Seyferts although their relative strengths may vary from galaxy to 
galaxy. 

Levenson et al. (2001) presented X-ray imaging and spectroscopy of a sample of 12 
Seyfert 2 galaxies and concluded that in order to fit the observed X-ray spectra it is required 
to combine the power-law Seyfert component with a thermal starburst emission. Jimenez- 
Bailon et al. (2005) presented XMM-Newton and Chandra observations of the Seyfert 2 
galaxy NGC 1808. They found the hard X-ray emission associated with an unresolved 
nuclear sources whereas the soft emission is dominated by a thermal component associated 
with an extended starburst. 

Terlevich et al. (1990) suggested to use the stellar Call triplet absorption feature in the 
IR continuum as a direct indicator on the presence of young unresolved stellar population 
in the nuclear regions of Seyfert galaxies. Heckman et al. (1997) and Gonzalez Delgado et 
al. (1998) found absorption line features associated with photospheres of O and B stars and 
their stellar winds in the ultraviolet and optical spectra of four Seyfert 2 galaxies: Mrk 477, 
NGC 7130, NGC 5135 and IC 3639 and thus presented direct evidence for the existence of 
nuclear starbursts in these galaxies. They found that the size of the nuclear starbursts in 
these galaxies ranges from several tens to a few hundred parsecs. The co-existence of W-R 
features in the optical and Ca II triplet in the near-IR part of the spectra, implies either a 
continuous star formation during more than ~ 10 Myr or two stellar generations with ages 
about 5-6 Myr and 10-20 Myr, respectively. Such starbursts are likely to drive the high- 
velocity outflows detected in the above and in Seyfert 2/starburst ultra luminous infrared 
galaxies (Gonzalez Delgado et al. 1998; Rupke et al. 2005). 

On the other hand, compact, bright stellar clusters or nuclear star clusters, were found 
in the centers of ~ 75% of local spirals and Virgo dwarf elliptical galaxies (Boker et al. 2002; 
Cote et al. 2006). Their radii (a few parsecs) are similar to those of globular clusters, however 
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they are 1-2 orders of magnitude brighter, more massive and may have complicated star 
formation histories with several episodes of star formation (Walcher et al. 2006). Ferrarese et 
al. (2006) claimed that massive galaxies (M gal > 10 10 M ) host supermassive BHs whereas 
less massive galaxies host only nuclear clusters. However Seth et al. (2007) presented 
evidences on the presence of super massive black holes in ~ 25% of galaxies which host 
nuclear star clusters. More than half of these galaxies (~ 15%) present a mixed AGN- 
starburst optical spectra and are classified as composite. 

Shlosman & Begelman (1989), Goodman (2003), Collin & Zahn (1999), Tan & Blackman 
(2005) have shown that accretion disks are gravitationally unstable outside of r ~ 10~ 2 pc 
and must fragment into self-gravitating clumps that eventually form stars. It was suggested 
that star formation reduces the rate of accretion and thus the luminosity of the central 
supermassive black hole (BH) by removing mass from the accretion flows and also due to 
radiative heating of the accretion discs. However non of these models took into consideration 
the negative feedback provided by the mechanical energy of the central starburst on the 
accretion flow. 

Thus circumnuclear star formation occurs at different space scales around the SMBH in 
many AGN galaxies. Here we note that the mechanical power of young nuclear starbursts 
might prevent through the cluster winds the accretion of interstellar matter from the bulges 
and disks of their host galaxies onto the central BHs. In such cases the BHs are fed with the 
matter injected by numerous stellar winds and SNe explosions that result from the multiple 
evolving sources. This implies that nuclear starbursts must strongly affect and perhaps even 
control the power of the central BH. In fact, it may be the dominant factor to be consider in 
order to understand the physics and relative contributions of the BH and starburst activity to 
the energy budget of the composite AGN/starburst galaxies. Lipari & Terlevich (2006) have 
incorporated different ingredients of this physics into their evolutionary unification scenario 
which seems to be able to explain many properties of AGNs and QSOs. 

The classic spherically-symmetric accretion model (Bondi, 1952; Frank et al. 2002) 
should then be modified if one is to apply it to the case of a massive BH at the center of a 
young stellar cluster. First, it should take into consideration the energy and mass supplied 
by massive stars within the star cluster volume. Second, it should account for radiative 
losses of energy from the hot thermal plasma. Third, the models should incorporate initial 
and boundary conditions as described by the star cluster wind theory (Chevalier & Clegg, 
1985; Canto et al. 2000; Silich et al. 2004; Tenorio-Tagle et al. 2007). 

Here we present a self-consistent semi-analytic theory of stationary spherically-symmetric 
flows driven by the thermalization of the kinetic energy supplied by massive stars inside mas- 
sive stellar clusters which includes the outflow of the injected matter and its accretion onto 
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a central massive BH. The solutions account for proper initial and boundary conditions for 
a variety of stellar clusters and black holes and the impact of strong radiative cooling on the 
dynamics of the thermalized injected matter. 

The paper is organized as follows. In section 2 we formulate our model and discuss the 
input physics and major simplifications. In section 3 we present a set of major equations. 
Boundary conditions and the selection of the proper solution from a family of integral curves 
are discussed in section 4. In section 5 we discuss the impact that a central BH provides on 
the star cluster driven flows. We discuss two hydrodynamic regimes separated in the star 
cluster mass vs star cluster radius diagram by a threshold line: the BH dominated regime 
below the threshold line and the radiative cooling dominated regime above the threshold line. 
In section 6 we use our model to calculate the accretion rates and BH accretion luminosities 
and compare them to those, predicted in the classical Bondi accretion theory. Section 7 
summarizes our results and gives our conclusions. 



2. The model 

Following Chevalier & Clegg (1985) we assume that massive stars are homogeneously 
distributed inside a spherical volume of radius Rsc and that the mechanical energy de- 
posited through stellar winds and supernovae, L S c, is thermalized via random collisions of 
the gaseous streams from neighboring sources. This results into a high temperature and a 
high thermal pressure that leads to a fast outflow of the injected matter, while composing a 
stationary star cluster wind. 

In presence of a massive, central black hole, a fraction of the deposited matter is to 
remain bound inside the cluster to eventually fall onto the center. This implies that in 
presence of a BH the stagnation point, the point where the expansion velocity, u w = km 
s -1 , is not at the star cluster center, as it happens if the cluster does not contain a BH and 
evolves in the quasi-adiabatic regime (see below), but instead is at a distance, R st , from the 
star cluster center. The position of the stagnation point thus becomes an important issue 
that defines both the upper limit for the accretion rate onto the central BH, and the amount 
of matter that the star cluster returns, in the form of a wind, to the interstellar medium. 

As it was mentioned by Nulsen & Fabian (2000), the dissipation of angular momentum 
is less of a problem for the hot plasma. Thus we restrict ourselves to spherically-symmetric 
solutions, despite realizing that some fraction of the injected material that remains bound 
inside the star cluster and that falls onto the center, eventually forms an accretion disk 
associated with the central black hole. We assume that this occurs at a sufficiently small 
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radii and neglect in this paper all effects associated with the redistribution of the residual 
angular momentum. 

We do not consider the relativistic effects and stop our integration at r = 3Rsh, the 
last stable orbit around the central black hole, where Rsh = 2GMbh/c 2 is the Schwarzschild 
radius, Mbh is the black hole mass, c is the speed of light and G is the gravitational constant. 

Throughout the calculations the parameters of the cluster, and the mass of the central 
black hole, are not allowed to change and thus only stationary solutions are discussed. 

Figure 1 presents a schematic representation of the internal structure of the bimodal 
flow induced by the presence of a central BH. In this case a fraction of the matter supplied 
by SNe and stellar winds is captured by the central BH and the rest is ejected as a wind into 
the ambient ISM. 



3. Main Equations 



The hydrodynamic equations for the steady-state, spherically symmetric flow that re- 
sults from the energy and mass deposition within young and massive star clusters with a 
radius Rsc, mass M$c and mechanical luminosity Lsc, are (see, for example, Johnson & 
Axford, 1971; Chevalier & Clegg, 1985; Canto et al. 2000; Silich et al. 2004): 
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for r < R S c and r > Rsc respectively. P, u, and p in equations (P - E]) are the thermal 
pressure, the velocity and the density of the thermalized matter. The mass and the energy 
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Fig. 1. — The structure of the flow that results from the thermalization of the supernova 
ejecta and stellar winds inside a young stellar cluster, when a massive black hole is located at 
the center. The radii of the internal and the external circles represent the stagnation radius 
R st and the star cluster radius R sc , respectively. The arrows indicate the direction of the 
flow. The black dot at the center marks the location of the black hole. 
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deposition rates per unit volume, q m = 3Msc /^Rsc anc ^ Qe = 3Lsc/^Rsci are assumed to 
be spatially constant inside the star cluster and equal to zero if r > Rsc- Q — n e n i-h(T, Z) 
is the cooling rate. A(T, Z) is the cooling function, T is the temperature and Z is the 
metallicity of the plasma. L sc = M sc V\ oa j r l, where Vaoo is the adiabatic outflow terminal 
speed. M(r) tracks the distribution of the stellar mass within the cluster, and Mbh is 
the mass of the central black hole. We have taken into consideration the gravitational pull 
provided by the central black hole and the star cluster and neglected the self-gravity of the 
reinserted gaseous component. 

One can easily integrate the mass conservation equation, both, inside and outside of the 
cluster, and rewrite equations ([T] - [6]) in the form: 
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in the region r > Rsc- c s = (7P/ p) 1 ^ 2 is the sound speed in the hot thermalized ejecta. The 
escape velocity, V esc , is: V esc = [2G(M(r) + M BH )/r] l l 2 if r < R sc , and V esc = [2G{M SC + 
Mbh) /r} 1 / 2 if r > Rsc, respectively. 

The presence of the BH does not affect the relation between the gas number density 
and the temperature at the stagnation point found in Silich et al. 2004: 
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is the adiabatic wind terminal speed, c s t and A(Z,T st ) are the 



sound speed and the cooling function calculated at r = R st . One can prove this result by 
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comparing the derivative of the expansion velocity at the stagnation point obtained from 
equation Qj with that obtained from equation ([7]), and requiring a finite derivative of density 
at the stagnation point. Note that Sarazin & White (1987) obtained similar relation from 
the energy conservation equation in their cooling flow model. 

In all calculations we use the equilibrium cooling function for optically thin plasma 
tabulated by Plewa (1995) and assume that the metallicity of the plasma is solar. 

4. Boundary conditions and the appropriate integral curve 

The thermalization of the mechanical energy supplied by massive stars within a young 
stellar cluster, causes a large thermal overpressure that drives away the injected matter in 
the form of a high velocity outflow - the star cluster wind. The smooth transition from 
a subsonic expansion of the high temperature thermalized ejecta, inside the star cluster 
volume, to the supersonic free wind outflow at r > Rsc, requires (see equations 191 and [T21) 
the sonic point (the point where the outflow velocity is equal to the local speed of sound) to 
be located at the star cluster surface (see Canto et al. 2000; Silich et al. 2004). Hereafter 
we will refer to this sonic point as the outer sonic point. 

In the case of stellar clusters with a central black hole, the gravitational pull of the 
BH prevents the escape of the injected matter from the central zones of the cluster and 
thus shifts the stagnation point from the star cluster center to a larger radius. In this 
case all mass continuously deposited by the cluster inside the central zone, limited by the 
stagnation radius, cannot escape from the gravitational well of the central BH and composes 
the accretion flow. The presence of the central BH results also in the existence of the second 
sonic point, between the stagnation radius and the star cluster center. To distinguish this 
sonic point from that at the star cluster surface we will refer to it as the inner sonic point. 

Thus the stagnation radius defines the upper limit to the accretion rate onto the central 
BH and also the fraction of mass that the cluster returns to the ambient ISM. This implies 
that the major problem that one has to solve in order to build a self-consistent hydrodynamic 
solution for the flow that results from the energy and mass deposition by a young stellar 
cluster with a central massive BH is reduced to the calculation of the stagnation radius. 

We show below that the proper position of the stagnation point is defined by the second 
boundary condition which is similar to that in the case of the Bondi accretion with 7 = 5/3. 
Specifically, we show that the inner sonic point must be located at the star cluster center. 
To avoid numerical problems associated with the central singularity we will assume that the 
inner sonic radius coincides instead with the last stable orbit associated with the central 
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black hole: 

Rsonic,in (1^) 

where R$h is the Schwarzschild radius of the central black hole. 

In order to select a proper integral curve we take a trial stagnation radius and then select 
T st from our first boundary condition which requires the outer sonic point to be located at 
the star cluster surface (see Silich et al. 2004; Tenorio-Tagle et al. 2007). We calculate then 
the number density of the plasma at the stagnation radius from equation ( fTBl and use these 
R s t,T st and n st as initial conditions to the backward integration from R st towards the star 
cluster center. 

Figure 2 presents the results of the integration for different values of the trial stagnation 
radius. If the selected R st is too large, the backward integration leads to a double- valued, 
unphysical solution, marked in Figure 2 by the dashed line. In this case the velocity turnoff 
point (marked by a cross symbol) coincides with the inner sonic point far from the cluster 
center. The turnoff point moves towards r = when the considered R st is smaller. This 
finally leads to an integral curve (marked in Figure 2 by the solid line) that approaches the 
last stable orbit at R = 3R s h with the sound speed. For even smaller R st the stagnation 
density is larger, the accretion velocity does not reach the sonic value at R — 3R s h and 
may even go to zero as it is shown in Figure 2 (dotted line). In order to learn how the inner 
boundary condition affects the solution, we have provided several runs with different values of 
R SO nic,in- Figure 3 shows that the stagnation radius is a weak function of R son ic,in and thus the 
inner boundary condition does not affect the solution significantly. We select as the proper 
solution the integral curve which has the second sonic point located at R son ic,in = 3Rsh- 

Figure 4 presents the distribution of the flow variables (velocity, number density and 
temperature) for a particular case of a 10 8 M black hole located at the center of a young 
stellar cluster whose mass M sc = 10 8 M and R sc = 40pc. In this case the stagnation 
radius, marked by the inner dotted line in panel a, is R s t = 2.7 pc. At larger radii the 
velocity grows almost linearly to reach the sonic value at the star cluster surface. It becomes 
supersonic outside the cluster and soon reaches the terminal value, V^, somewhat smaller 
than the adiabatic terminal speed value, as radiative losses deplete some energy inside the 
cluster and in the free wind region. 

In the region between the stagnation radius and the black hole, the matter deposited 
by stellar winds and supernovae composes a stationary accretion flow. The absolute value 
of the velocity grows rapidly in this region. However the flow remains subsonic as radiative 
losses are not able to compensate the heating of the in-flowing plasma. This leads to a rapid 
increase of temperature (panel c) despite the increase in density and thus of cooling of the 
in-falling matter, as it approaches the star cluster center (panel b). 
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log (r/Rs h ) 

Fig. 2. — Possible integral curves. Three possible integral curves marked by dashed, solid and 
dotted lines correspond to different stagnation radii: 5 pc, 3.3 pc and f pc, respectively. The 
calculations assumed a star cluster mass M sc = 1O 8 M , a radius, R S c = 100 pc, a black 
hole mass, M B h = 10 8 M Q , an adiabatic wind terminal speed, Va j00 — 1500 km s" 1 , and 
solar metallicity. The dashed line presents the unphysical double-valued solution. The solid 
line shows the selected solution which satisfies for both boundary conditions. The dotted 
line presents another unphysical branch of integral curves which tends towards positive flow 
velocities around the black hole. The normalization velocities are uq = 10 4 km s" 1 , for 
dashed and solid lines, and u = 10 2 km s~ x for the dotted line, respectively. 
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Fig. 3. — The impact of inner boundary condition on the solution. The cal- 
culations were provided for five different values of inner sonic radius: R SO nic,in = 
3R S h', 10.5Rsh', 98R S h', 1009 Rsh and 2385R S h- The star cluster and black hole parameters 
are identical to those in Figure 2: M S c = 1O 8 M , R S c = 100 pc, Va j00 — 1500 km s _1 and 
M B h = 1O 8 M . Cross symbols represent the results of the calculations. 
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Radius (pc) 

Fig. 4. — The structure of the flow in the case of a young stellar cluster with a central massive 
black hole. Panels a, b and c present the distribution of the flow velocity, number density 
and temperature, respectively. The calculations assumed a 10 8 M black hole located at 
the center of a young stellar cluster whose mass is M sc = 10 8 M and R sc = 40pc. The 
adopted mechanical luminosity of the cluster is 3 x 10 42 erg s -1 . This value corresponds to 
the average mechanical luminosity of a young stellar cluster with a Salpeter initial mass 
function (Leitherer et al. 1999). It was assumed that the adiabatic wind terminal speed is, 
Va,oo = 1500 km s -1 . Vertical dotted lines in panel a mark the stagnation radius, R st , and 
the star cluster radius, Rsc, respectively. 
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Figure 5 shows how the location of stagnation point depends on the mass of the black 
hole and that of the stellar cluster, when the radius of the cluster remains fixed. R st be- 
comes larger when the gravitational pull becomes stronger either because one considers more 
massive black holes or clusters. 
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Fig. 5. — The position of the stagnation point. Panel a shows how the stagnation radius 
grows with the mass of the black hole for a fixed (10 8 M Q ) star cluster mass, whereas panel 
b shows the value of the stagnation radius as a function of the star cluster mass for a fixed 
(1O 8 M ) black hole mass. In all cases the radius of the cluster is 40 pc. 



5. The threshold mechanical luminosity 

Silich et al. (2004), Tenorio-Tagle et al. (2007) and Wiinsch et al. (2008) have thor- 
oughly discussed the impact of radiative cooling on the inner structure of star cluster driven 
flows in the case without BH and without accounting for the contribution to the gravitational 
field imposed by the cluster. They showed that strong radiative cooling changes drastically 
the pressure gradient in the inner zones of compact and massive star clusters and found a 
critical, threshold mechanical luminosity which separates, in the L$c - Rsc parameter space, 
clusters with R st = pc from those evolving in the bimodal, catastrophic cooling regime. 
In the latter regime, strong radiative cooling promotes the displacement of the stagnation 
point out of the star cluster center and leads to the accumulation of the matter injected 
inside the stagnation volume as it becomes thermally unstable, while the outer zones of the 
cluster drive a stationary outflow. On the other hand, as shown above, stellar clusters with a 
central BH evolve always in a bimodal regime because the gravitational force goes to infinity 
when the distance to the BH goes to zero. This establishes the position of a stagnation point 
within the flow and inhibits the escape of the hot plasma from the cluster inner zones. The 



implication is thus that the two physical processes here considered: the gravitational pull 
from the central massive BH and strong radiative cooling, both promote the existence of a 
stagnation point within the flow and thus both compete in defining the final hydrodynamic 
solution. 

In order to calculate the threshold energy, we choose the stagnation temperature which 
leads to the maximum thermal pressure at the stagnation point. This is defined by the 
condition that dP st /dT st = which, together with equation (TT3|) . yields 



boundary conditions are fulfilled, e.g. until the outer sonic point accommodates at the star 
cluster surface and the inner one reaches the 3Rsh radius. 

Figure 6 displays the threshold mechanical luminosity calculated for stellar clusters of 
different radii which contain black holes of different masses at their centers. The gravitational 
pull of the BH does not affect too much the value of the threshold luminosity in the case of 
large clusters. Although in all of these cases the critical luminosity is several times smaller 
than in the case without a BH (see Figure 6). 

The gravitational pull of the BH becomes progressively more important for more com- 
pact clusters. The threshold lines turn up when the considered BH mass becomes comparable 
to the critical cluster mass (see Figure 6). In the case of even more compact clusters, cooling 
cannot compete with gravity. The position of the stagnation point is defined then by the 
mass of the BH and the stationary solution exists regardless of the cluster mass. Perhaps 
the only limitation here arises when R st becomes larger than Rsc- This implies that the 
critical luminosity does not exist below Rsc,crit- The critical radii, Rsc,crit, depend on the 
mass of the BH and are marked in Figure 6 by thin vertical lines for different values of the 
BH mass. 

If the cluster does not contain a BH and its gravitational pull is negligible, the stagnation 
radius is defined by the excess mechanical luminosity, L$c, over the threshold value (Wiinsch 
et al. 2007): 



From this it is clear that the impact of radiative cooling on the inner structure of 
the flow is to become progressively more important for star clusters with larger mechanical 
luminosities (proportional to the star cluster mass). For low mass clusters, the gravitational 
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Then we iterate the stagnation radius and the star cluster mechanical luminosity until both 
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Fig. 6. — The threshold mechanical luminosity. The thick solid line marks the threshold 
luminosity for cases without a central black hole. The dotted, dash-dotted, dashed and thin 
solid lines display the threshold luminosity for clusters with 10 7 M 10 8 M Q 5 x 10 8 M and 
10 9 M Q black holes, respectively. Thin vertical lines mark critical radii, R,sc,crit- Cooling 
cannot compete with gravity if the star cluster radius, Rsc < Rsc,crit- In this parameter 
space gravity defines the stagnation radius regardless of the star cluster mass and the critical 
luminosity does not exist. It was assumed an adiabatic wind terminal speed equal to 1500 km 
s _1 and solar metallicity. 
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pull of the black hole is to dominate and define the position of the stagnation point, but 
radiative cooling is to become an important factor when one considers more massive clusters. 

We expect that above the threshold line the internal structure of the flow would present 
three distinct zones, as displayed in Figure 7. The outer zone would conform a quasi- 
stationary outflow. This runs from an outer stagnation radius R s t,cooi, set by strong radiative 
cooling, to meet the sound speed at the cluster surface. The central region will enclose the 
accretion flow onto the BH. The outer boundary of this flow would be the stagnation radius 
defined by the BH (R s t,BH)- Between these two stationary flows the thermalized ejecta is 
thermally unstable and rapidly decays into two phases: a hot plasma whose parameters are 
similar to those found at the outer stagnation point, and a collection dense clouds which 
result from the thermally unstable plasma and are completely or partially photo-ionized by 
the stellar UV and BH hard radiation (Wiinsch et al. 2008). These clouds can fall onto 
the BH or become gravitationally unstable and support some level of star formation inside 
such a cluster. Numerical calculations and a thorough discussion of this regime will be the 
subject of a forthcoming communication. 



6. Accretion rates and BH luminosities 

For clusters below the threshold line, the mass deposited by stellar winds and supernovae 
explosions inside R st is not able to escape the stagnation volume and in the stationary regime 
has to fall to the center providing fuel to the central BH. Thus in the stationary regime the 
stagnation radius, R st , defines the accretion rate onto the central object: 

M acc = M sc (^j. (17) 

In this respect it is instructive to compare our theory with the classic spherically- 
symmetric accretion theory (Bondi, 1952; see Frank et al. 2002 and references therein). 
If the polytropic index 7 = 5/3, the Bondi accretion rate is (Frank et al. 2002): 

M B = K&M* H q™, (18) 

C ISM 

where Pism and cism are the density of the ISM and the speed of sound at infinity, respec- 
tively. We associate these quantities with Chevalier & Clegg's central values (see Canto et 
al. 2000): 



(19) 
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Fig. 7. — The expected velocity pattern of the flow in the case of clusters located above the 
threshold line. The stationary solution does not exist in this case. Positions of the inner and 
outer stagnation points are defined by the gravitational field of the central BH and strong 
radiative cooling, respectively. The intermediate zone is thermally unstable. 
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Figure 8 compares our semi-analytic results (equation [TTj) with Bondi accretion rates 
(equation [T8j) . The calculations were performed for a set of clusters which contain a 10 s M 
BH at the center, all clusters have the same radius of Rsc = 40 pc. It was assumed that the 
adiabatic wind terminal speed is Vaoo = 1500 km s -1 and a plasma with solar metallicity. 
In order to relate the mechanical luminosity of the stellar cluster to the corresponding star 
cluster mass we have used a relation which approximates the results of the Starburst 99 
synthesis model for young stellar clusters (Leitherer et al. 1999): 

3 x 1040 (i^) < 22 > 

The star cluster mechanical luminosity has been normalized to the critical value derived 
in the previous section. Thus when Lsc/L crit <C 1, clusters evolves in the quasi-adiabatic 
regime, whereas when L sc lL crit — > 1 the star cluster parameters approach the threshold 
values. The semi-analytic results are marked by the cross symbols and the accretion rates 
predicted by the Bondi equation are shown by the solid line. 

Figure 8 shows that well below the threshold line (when Lsc/L cr u < 0.1) Bondi's 
formula (equation [T8|) is in good agreement with our numeric results. For such clusters one 
can approximate the spherically-symmetric accretion rate onto a BH located at the center 
of a young massive cluster with Bondi's expression: 

M acc = nG 2 M 2 BH ^, (23) 

where p c and c c are taken from equations ffl~9l) - fl2Tl) . One can then use equation (|23|) in order 
to obtain an analytic expression for the stagnation radius in this (Lsc ^ L cr i t ) parameter 
space. Indeed, the mass accretion rate is 

Macc = 4nq m R 3 st /3, (24) 

where the mass deposition rate per unit volume, q m , is 
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Fig. 8. — The comparison of the approximate analytic formulae with numeric accretion rates. 
Cross symbols present the numeric results for a 10 8 M BH embedded into stellar clusters 
of different masses, all having radius of Rsc — 40 pc. Solid line displays the accretion rates 
calculated in the Bondi approximation with the adiabatic wind input parameters. 
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Combining equations (1241) and fl25l) with equation (1231) . one can obtain: 

RM ~ { 2L SC 4 ) Rsc - (26) 

For star clusters whose mechanical luminosities, Lsc, are comparable to the threshold 
value, L crit , the calculated accretion rates exceed substantially those predicted by Bondi's 
equation. In this case one has to use the semi-analytic model in order to find the size of the 
stagnation zone and then equation (ITT)) in order to derive the accretion rate. 

Having the accretion rates and adopting the normal accretion efficiency, 77 = 0.1, one 
can calculate the BH luminosity. Figure 9 compares the calculated BH luminosities with the 
Eddington limit, L Edd = 1.3 x 10 38 M BH /M Q erg s -1 , for star clusters evolving in different 
hydrodynamic regimes. Lbh grows in the case of more compact clusters. However, it remains 
well below the Eddington limit even when the star cluster mechanical luminosity reaches the 
threshold value (Figure 9, panel a). The BH luminosity can approach the Eddington limit 
(see Figure 9, panel b) in the case of very compact clusters whose radii are smaller than the 
critical values, Rsc,crit, marked in Figure 6 by vertical lines. This puts a limit for our model 
in the case of very compact clusters. 




Fig. 9. — The BH accretion luminosity. Panel a presents the accretion luminosity of a 
10 s M Q BH embedded into different star clusters. Solid, dashes and dotted lines display the 
results of the calculations for star clusters whose radii are Rsc = 40 pc, Rsc = 30 pc and 
Rsc = 10 pc, respectively. The last point on every line displays the BH accretion luminosity 
when the star cluster reaches the threshold line. Panel b presents the accretion luminosity 
of a 10 8 M BH embedded into very compact (Rsc — 3 pc) star clusters, all located on the 
vertical dot-dashed line in Figure 6. 
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7. Conclusions 

We have developed a self-consistent, stationary solution for spherically symmetric ac- 
cretion flows which are formed inside young stellar clusters with a central supermassive BH. 

We have shown that the thermalization of the kinetic energy released by massive stars 
inside young stellar cluster results in a bimodal solution which presents an accretion of the 
injected matter onto a central BH in the inner zones of the cluster, and the ejection of the 
deposited matter from the outer zones of the cluster in the form of the fast superwind. We 
suggest that superwinds prevent the accretion of the ambient interstellar gas from the bulges 
and disks of their host galaxies onto the central BHs and that in such cases the BHs are 
fed with the matter re-inserted by massive stars in the form of numerous stellar winds and 
SNe explosions. The accretion rate and the BH luminosity are then defined by the central 
starburst but not by the gravity and the interstellar gas distribution in a host galaxy. 

The hydrodynamics of the accreted and ejected matter depend on the location of the 
stellar cluster hosting a BH in the cluster mass - black hole mass - star cluster radius 
parameter space. There is a surface in this 3D parameter space which separates clusters 
evolving in the stationary regime from those which cannot fulfill stationary conditions. If 
the mechanical luminosity of the cluster exceeds the critical value, a hot plasma inside the 
cluster is thermally unstable. The flow is highly non-stationary and presents a complicate 
velocity patten with an outer stagnation point defined by strong radiative cooling and the 
inner one whose position is defined by the mass of the BH. The hydrodynamical structure 
and time evolution of the resulting flow in such cases must be calculated numerically. 

Clusters whose mechanical luminosity is smaller than the critical value, compose sta- 
tionary accretion flow in the central zones and form stationary outflows, the star cluster 
winds, in the outer zones of the cluster. 

We used our model to calculate the accretion rates and the accretion luminosities of 
BHs at the centeres of young star forming regions. The classic, Bondi's accretion theory 
shows a good agreement with our semi-analytic model, but only in the case of low mass 
clusters located well below (L$c < 0.1L crit ) the critical luminosity in the L$c — Mbh — Rsc 
parameter space. Thus one has to use a semi-analytic approach in order to calculate the 
accretion rates and BH luminosities in the case of more energetic clusters. 

In the case of extended starbursts, the BH luminosities fall well below the Eddington 
limit. However, the accretion luminosities grow rapidly for more compact clusters and for 
very compact clusters can approach the Eddington luminosity, as shown in Figure 9, panel 
b. 
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The model, here developed, is required in order to advance our knowledge regard- 
ing the relative contributions of supermassive BHs and central star bursts in composite, 
AGN/starburst galaxies and will be used for interpretation of observational properties of 
such objects in a further communication. 

We thank our anonymous referee for valuable comments and suggestions. SS wishes 
to express his thanks to Isaac Shlosman for many useful discussions during his stay in the 
University of Kentucky where this study was designed. This project has been supported by 
CONACYT - Mexico research grants 47534-F and 60333. 
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